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Abstract 

Large-scale storage is a promising emerging technology to realize a reliable smart-grid since it 
can enhance sustainability, reliability and asset utilization. At fast time-scales, storage can be used to 
mitigate fluctuations and time-variations on renewable generation and demand. On slower time scales, 
it can be used for load shifting and reducing the generation costs. This paper deals with the latter case 
and studies the following important decision problem; Given a network and an available storage budget, 
how should we optimally place, size and control the energy storage units across the network. When the 
generation cost is a nondecreasing convex function, our main result states that it is it is always optimal 
to allocate zero storage capacity at generator buses that connect to the rest of the power grid via single 
links, regardless of demand profiles and other network parameters. This provides a sound investment 
strategy for network planners, especially for distribution systems and isolated transmission networks. 
For specific network topologies, we also characterize the dependence of the optimal production cost on 
the available storage resources, generation capacities and flow constraints. 

I. Introduction 

Emerging energy storage technologies in the electricity grid can be potentially used for various 
services, e.g., to enhance sustainability, reliability, efficiency and asset utilization [[Il-[|4l. The 
economic value [I51-0 and national impact [HI, flU of these services has been emphasized in 
various recent studies. Over the last decade or so, these technologies have shown remarkable 
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technical improvements and their cost of production has dropped significantly lfT0l - [fT2l . Also, 
the obvious challenges in integrating intermittent renewable generation has spurred the interest 
in large-scale storage to deal with their variability [ir3l - [fT6ll . In summary, the authors in [[8] 
identify the use of energy storage as "critical to achieving national energy policy objectives 
and creating a modern and secure electric grid system." For this promising new technology, we 
study an efficient investment strategy of placing and sizing such units in a network to minimize 
generation cost. 

Applications of energy storage can be classified according to time scale of operation 0, [|5l, 
ifTTl . At fast-time scales (seconds to minutes), storage is mainly used to mitigate variability of 
renewable generation and demand and reduce the role of costly ancillary services to balance 
demand and supply, e.g., lfT8l - ll2T1l . On a slower time scale (over hours), bulk storage devices 
aim at load shifting, i.e., generate when it is cheap and supply the residual load using storage 
dynamics ifTOl . In the present work, we only deal with the latter. Our results are therefore 
complementary to those of lITSl - llITI . The short-time scale strategies developed there can be 
implemented " on top " of the long-term strategies developed here. 

Now, we provide a brief overview of the relevant literature. Optimal control policy for storage 
units has been extensively studied. While the authors in [|20ll . [|22ll examine the control of a single 
storage device without a network, [1211 . [|23l . [|24l explicitly model the role of the networks in 
the operation of distributed storage resources. The engineering constraints are designed based 
on one of the two popular power flow models [|25l - [|27l . namely, AC power flow and DC power 
flow. Storage resources at each node in the network are assumed to be known a priori in these 
settings. 

For the emerging smart grid technology, the investment decision problem of selecting, sizing 
and placement of distributed energy resources is gaining importance. As noted in [|28l . it is 
critical to "assess the technical and economic attributes of energy storage specifically reflecting 
the operational demands and opportunities presented by the smart grid environment. Without the 
ability to analyze network features of transmission and distribution, storage system technologies 
and their efficiencies, along with their cost beneflts for various value streams, there is no ability 
for the utility to make comparative business decisions that will enable the optimal siting of 
energy storage." Several researchers have analyzed aspects of this problem, e.g., ||29l , [[30l 
using purely economic arguments, without explicitly considering the network constraints of the 



physical system. Authors in [|20l . [|3T| have looked at optimal sizing of storage devices in single- 
bus power systems, while Kanoria et al. ||2T1| compute the effect of sizing of distributed storage 
resources on generation cost for specific networks. Gayme et. al [|23l study a similar problem in 
IEEE benchmark systems [[32l . Recently, a more general framework to study the optimal storage 
placement problem in generic networks has been formulated and studied through simulations in 
ll33l . [|34l using tools similar to optimal power flow (OPE) . The resulting problem is non-convex 
and hence Sjodin et al. in ll33ll use the approximate DC-OPE [|26l . Il27l while Bose et al. in ||34| 
use the relaxation of AC-OPE [|35l - [l38l based on semidefinite programming [|39l , [|40l . 

In this paper, we consider the investment decision problem of how much storage to place 
on each node in a network given an available storage budget. As a by-product, we also derive 
the optimal control policy for the storage units. Our main result shows that when minimizing 
a convex and nondecreasing generation cost with any fixed available storage budget, there 
always exists an optimal storage allocation that assigns zero storage at those generator nodes 
that connect via single lines to the network, for arbitrary demand profiles and other network 
parameters. It suggests that in most distribution networks and isolated transmission networks, it 
is always optimal to allocate the entire available storage budget among demand buses. Also, this 
work provides an analytic justification to the conjecture that the optimal allocation of storage 
resources mainly depends on the network structure as opposed to the total available storage 
budget as noted through simulations in [|23l . [|33l . [|34l . Eor generator buses with multiple line 
connections, however, we show that the optimal storage placement may not in general place zero 
storage capacity at such nodes. Eurthermore, we study the dependence of optimal production 
cost on the available storage budget, line flow and generation capacities for specific network 
topologies. The physical network has been modeled using DC power flow. Preliminary versions 
of our result have appeared in PTI and 02). 

The paper is organized as follows. We formulate the optimal storage placement problem in 
Section |Ill The main result is proven and discussed in Section Hill More detailed analysis of 
the problem on specific network topologies is discussed in Section |IVl Concluding remarks and 
directions for future work are outlined in Section Technical details appear in the Appendix. 



II. Problem formulation 

Consider a power network that is defined by an undirected connected graph Q on n nodes (or 
buses) J\f = {1,2, For two nodes k and / in J\f, let A; ~ / denote that k is connected to / 
in ^ by a transmission line. For the power flow model, we use the linearized DC approximation; 
see [|26l for a detailed survey. In this approximation, the network is assumed to be lossless, the 
voltage magnitudes are assumed to be at their nominal values at all buses and the voltage phase 
angle differences are small. 

Time is discrete and is indexed by t. Now consider the following notation. 

• dk{t) is the known real power demand at bus k E M at time t. Demand profiles often show 
diurnal variations [|43l . i.e., they exhibit cyclic behavior. Let T time-steps denote the cycle 
length of the variation. In particular, for all k E J\f, t > 0, assume 

4(t + T) = 4W. 

• gk{t) is the real power generation at bus k E J\f at time t and it satisfies 

0<9k{t)<g,, (1) 

where, gg. is the generation capacity at bus k. 

• Cfc (gk) denotes the cost of generating power gk at bus k E M. The cost of generation is 
assumed to be independent of time t and depends only on the generation technology at 
bus k. Also, suppose that the function Ck : ffi^ — t- is non-decreasing, differentiable and 
convex. A cost function commonly found in the literature [|23l . [|25l . (371, ll38l that satisfies 
the above assumptions is: 

Cfc {gk) = ik,2gl + ik,igk + 7fc,o, 

where 7fc,2, 7fc,i, 7fc,o are known nonnegative coefficients. 

• 9k{t) is the voltage phase angle at bus k E M at time t. 

• For two nodes k I, let pki{t) be the power flowing from bus k to bus / at time t. It 
satisfies 

Pki{t) = Vki [Ok{t) - ei{t)] , (2) 



where i/ki is the admittance of the Une joining buses k and /. Also, the power delivered 
over this line is limited by thermal effects and stability constraints and hence, we have 



\Pki{t)\ < fkh 



(3) 



where f^i is the capacity of the corresponding line. 
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Fig. 1: Power balance at node k E M. 



• 7fc(t) and 5k{t) are the average charging and discharging powers of the storage unit at bus 
k E M dX time t, respectively. The energy transacted over a time-step is converted to power 
units by dividing it by the length of the time-step. This transformation conveniently allows 
us to formulate the problem in units of power [|34l . Let < a^, < 1 denote the charging 
and discharging efficiencies, respectively of the storage technology used, i.e., the power 
flowing in and out of the storage device at node k E M aX time t is a^7fc(t) and -^^5k{t), 
respectively [|20l . [|44l . The roundtrip efficiency of this storage technology is a = a^as < 1. 

• Sfc(t) denotes the storage level at node k E J\f at time t and s° is the storage level at node 
k at time t = 0. From the definitions above, we have that 

Sk{t) = s° + 5^ (^a^lkir) - ^4(r)^ . (4) 

For each k E J\f, assume s° = 0, so that the storage units are empty at installation time. 



6fc > is the storage capacity at bus k. Thus, Skif) for all t satisfies the following: 



< Skit) < hk. (5) 

• h h the available storage budget and denotes the total amount of storage capacity that can 
be installed in the network. Our optimization algorithm decides the allocation of storage 
capacity at each node k e Af and thus, we have 

Y,h<h. (6) 

• The charging and discharging rates of each storage device are assumed to be upper-bounded 
by ramp limits. Suppose these limits are proportional to the capacity of the corresponding 
device, i.e., for all k E M, 

< -ik{t) < e^bk, (7a) 
< Skit) < esh, (7b) 

where G (0, ^] and es G (0, as] are fixed constants. 
Balancing power that flows in and out of bus k E M at time t, as shown in Figure [B we have: 

9kit) - dkit) - -fkit) + Sk it) = J2 Pki (t) ■ (8) 

l~k 

Also, optimally placing storage over an infinite horizon is equivalent to solving this problem 
over a singe cycle, provided the state of the storage levels at the end of a cycle is the same as 
its initial condition ||34|| . Thus, for each k G A/", we have 

(a77fcW - =0. (9) 

For convenience, denote [T] := {1, 2, . . . , T}. Using the above notation, we define the following 
optimization problem. 



Storage placement problem P: 

T 

minimize ^^^0^(5^^ {t) ) 

k£Af t=l 

over {gk (t) , 7fc (t) , 4 (t) ,9 kit), p^i [t] ,bk), 

keAf, kr^l, te [T], 

subject to O,©,©,®,®,®,©,®,®, 

where, ([1]) represents generation constraints, (|2),(|3) represent power flow constraints, (IS]),©,©,® 
represent the constraints imposed on the charging/discharging control policy of the energy storage 
devices, ® represents the power balance constraints at each bus of the network and ® represents 
the constraint on the sum of the capacities of all storage devices being no greater than the 
available storage budget. 

Restrict attention to network topologies where each bus either has generation or load but 
not both. Partition the set of buses J\f into two groups Ac and A//? where they represent the 
generation-only and load-only buses respectively and assume A/g and Md are non-empty. For 
any subset JC of Ag, define the following optimization problem. 
Restricted storage placement problem n'^: 

T 

minimize '^'^^k {gk (t) ) 

k&N t=l 

over {gk (t) , 7fc (t) , 4 (t) ,Bk{t), pki (t) ,bk), 

keAf, kr^l, te [T], 
subject to 0,©,©,@,®,®,©,®,®, 
bi = 0, i eK,. 

Problem II'^ corresponds to placing no storage at the (generation) buses of the network in subset 
/C. We study the relation between the problems P and II'^ in the rest of the paper. 

We say bus k e Af has a single connection if it has exactly one neighboring node / ~ k. 
Similarly, a bus k e Af has multiple connections if it has more than one neighboring node 
in Q. We illustrate the notation using the network in Figure |2] Afc = {1,2,7} and Afo = 
{3,4,5,6}. Buses 1 and 2 have single connections and all other buses in the network have 



multiple connections. 
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Fig. 2: A sample network. 
III. Main Result 

For a subset /C C Ag, let and txJ^ be the optimal values for problems P and II'^, 
respectively. Now, we are ready to present the main result of this paper. 

Theorem 1. Suppose K, C Afc and each node i E K. has a single connection. If P is feasible, 
then n'^ is feasible and = nf. 

The theorem states that, for any available storage budget, there always exists an optimal 
allocation of storage capacities that places no storage at any subset of generation buses with single 
connections. The result holds, regardless of the demand profiles and other network parameters, 
such as line flow constraints and impedances of transmission lines. In our model, the demand 
profiles of the load buses are deterministic. However, since Theorem [H holds for arbitrary demand 
profiles, our result extends to the case where the load is stochastic. 

Note that we restrict attention to generator buses in JC that have single connections only. This 
is applicable in many practical scenarios, as discussed in Section IIII-B[ The result is not true, 
in general, if /C includes generator buses with multiple connections. We provide an example in 
Section IIII-CI to illustrate this. 

A. Proof of Theorem [7] 

We only prove for the case where the round-trip efficiency is a < 1, but the result holds for 
a = 1 as well. Assume P is feasible throughout. For any variable z in problem P, let z* be the 
value of the corresponding variable at the optimum. In our proof, we use the following technical 
result. 



Lemma 1. Suppose : R — ?> R convex and dijferentiable. Then, for any xi < and 

< r] < {x2 — Xi): 

(f){Xi +7]) + (f){x2 - < (piXl) + 0(X2). 

Proof: We only prove the result for < < (x2 — xi)/2. The other cases can be shown 
similarly. The function 0(-) is continuous and differentiable. From mean-value Theorem [[39ll . it 
follows that there exists G (xi, xi + rf) and .^2 G (x2 — ?7, X2), such that + 77) — 0(xi)] — 
[(t){x2) - (f){x2 - T])] = 4>'{^i)r] - 4>'{^2)v- Also, ^1 < ^2 and hence 4>'{^i) < since is 

convex. This completes the proof. ■ 
Consider node z G /C and j ^ i. Node j is uniquely defined as i has a single connection. 
It can be shown that problem P, in general, has multiple optima. In the following result, we 
characterize only a subset of these optima. 

Lemma 2. There exists an optimal solution of P such that for all t G [T] and all i ^ IC, j ^ i, 

(a) g*{t)j:mit) = 0, 

(b) gm < 

The first part of Lemma |2] essentially says that for some optimum solution of P, the storage 
units should not charge and discharge at the same time step if there is positive generation at 
the same bus at that time step. This is expected since the round-trip efficiency of the storage 
devices a = a^as is less than one and since the generation cost is a nondecreasing function. 
The second part can be interpreted as follows. Power that flows from bus i to bus j at each 
t G [T] is Pij{t) = gi{t) — '^i{t) + 5i{t) and we have Pij{t) < fij. But Lemma lIHb]) states that 
there exists an optimum for which, g*{t),t G [T] itself defines a feasible flow over this line. 

Proof: The feasible set of problem P is a bounded Q polytope and the objective function 
is a continuous convex function. Hence the set of the optima of P is a convex and compact 
set ll39l . Now, with every point in the set of optimal solutions of P, consider the function 
^jg;cig[T] ilii't) + ^iii))- This is a linear and hence continuous function on the compact set of 
optima of P and hence attains a minimum. Consider the optimum of P where this minimum is 



' Without loss of generality, let bus 1 be the slack bus and hence 6i(t) — for all f G [T]. Boundedness of the set of feasible 
solutions of P then follows from the relations in O, (|2j, (H}, (|6) and (|7},. 



attained. We show that for this optimum, g*it)-f*{t)6*{t) = and g*{t) < fij for all t G [T] and 
z G /C, j ~ z. 

(a) Suppose, on the contrary, we have 5'*(to) > 0, 7* (to) > and 5* (to) > for some to ^ [T]. 
Define 

Ag' := min |(1 - «)7;(to) , ^^K^o) ,^?;(to)| • 

Note that Ag' > 0. Now, for bus i, construct modified generation, charging and discharging 
profiles gi{t),6i{t),%{t),t G [T] that differ from g* (t) , 6* (t) , ^y* (t) only at to as follows: 

^^(to) ■.= g:{to)-Ag', 

%ito) ■.= i:ito)--^Ag', 
1 — a 

Uto) := SKto) - j^Ag'. 

1 — a 

Note that, for all t G [T], the storage level Si{t) and the power Pij{t) flowing from bus i 
to bus j remain unchanged throughout. It can be checked that the modified profiles define 
a feasible point of P. Since Cj(-) is non-decreasing, we have Q (^j(to)) < Ci{g*{to)) and 
hence the additivity of the objective in P over i and t implies that this feasible point has an 
objective function value of at most p^. It follows that this feasible point defines an optimal 
point of P. However, we have 7j(to) + 5i(to) < 7j*(to) + <^i'(to) and thus, this optimum of P 
has a strictly lower J2ieict£[T] (7i(^) + '^i(^))' contradicting our hypothesis. This completes 
the proof of g*{toh:{to)6*{to) = 0. 

(b) If g*{t) = for all t G [T], then g*{t) < fij clearly holds. Henceforth, assume maxtgpr] g*{t) > 
0, and consider any to G [T], such that 5'*(to) = m.axt^[x] fl'*(t)- 

If 7* (to) = 0, then, 

max^*(t) = ^*(to) 

t€lT] 



=pUto)+i*{to)-s:{to) 

< (10) 



and Lemma 12© holds. 



Suppose now that 7* (to) > and hence 5* (to) = from Lemma |2l|al). First, we show that 
the storage device discharges at some point after tg. 



>o >o 

We also have s*(T) = s° = by hypothesis. Thus the storage device at node i needs to 
discharge in [to + 1, T] and hence a^j*{t) - ^5,*(t) < for some t G [to + 1, T]. Let ti be 
the first time instant after to when the storage device at bus i is discharged, i.e. 

ti :=min jtG [to + l,T] | a^7*(t) - — 5*(t) < ol. (11) 
t as ) 

Thus, 5*{ti) > 0. Define 

:= min |7;(to) , ^5*(ti) , (7*(to)|. (12) 

Then A^f > 0. Now, consider the case where: 

^;(ti)>0, and g:{to)<g:{ti) + aAg. (13) 

Since g*{ti) > and 6*{ti) > 0, then 7*(ti) = 0, by Lemma [2llaD. In that case, g*(ti) + 
6*{ti) = p*j{ti) is the power that flows from bus i to bus j at time ti. Combining (fT2l) and 
(fT3l) . we have 

max^*(t) = ^*(to) 
£G[r] 

<g:{h) + aAg 
<g*{h) + 6*{h) 

Hence, Lemma [2110) holds when (fT3l) is satisfied. Next, we show that if (fT3l) does not hold, 
then we can construct an optimum of P with a lower J2i£Kt&[T] (7«(^) + <^i(t)) and this 
contradicts our hypothesis. 

Suppose (fT3l) does not hold. If g*{ti) = 0, then we have 



g*{to) >Ag> aAg = g*{h) + aAg. 



Thus, it suffices to only consider tlie following case: 



gKt^)>gKh) + a^g. (14) 

Construct the modified generation, charging and discharging profiles at node i, gi{t), 5i{t),'^i{t) 
using (fT2l) . that differ from g* (t) , 6* (t) , •y* (t) only at to and ti as follows: 

gM = gKto) - Ag, Uh) = g*{h) + a^g. 

l^{to) = y*{to)-Ag, 7,(ti)=7;(tO, 
Uto) = S*{to) = 0, Uti) = - aAg. 

Also, define the modified storage level Sj(t) using 7j(t) and Si{t). To provide intuition to the 
above modification, we essentially generate and store less at time to by an amount Ag. This 
means at a future time ti, we can discharge aAg less from the storage device and hence 
have to generate aAg more to compensate. To check feasibility, it follows from (fT2l) . that 

for t = to,ti, we have 

0<g^{t)<g^, 

< m < ^,bi 

< Ut) < ^sb*. 

Also, the line flows Pij{t) remain unchanged. For the storage levels, it can be checked that 
the following holds: 

< s*{to - 1) < < s*it) < b*, for t e [to,ti - 1], 

Si{t) = s*(t), otherwise. 

This proves that the modified profiles define a feasible point for P. The cost satisfies 



Ci {gt{to)) + Ci (giih)) 

< Ci {g*{to) - aAg) + q {g*{ti) + aAg) 

<cdgKto)) + cagKh)). 



(15a) 
(15b) 



Equation (II Sal ) follows from the non-decreasing nature of Cj(-) and equation (|15b|) follows 
from using ([14] ) and Lemma [T] Thus the modified profiles gi{t) , 5i{t) , ^i{t) define a feasible 
point of P with a cost at most j>^, and, hence, are optimal for P. However, we also have 

>0 

Thus, the modified profiles define an optimum of P with a lower Ylii^K.t<^[T] (7«(^) + ^j(^))- 
This is a contradiction and completes the proof of the Lemma. 

■ 

To prove Theorem [1] consider the optimal solution of P that satisfies Lemma [2©. For all 
i E IC, g*{t) itself defines a feasible flow over the line joining buses i and j, where j is the unique 
neighboring node of i. Now the proof idea is as follows. For i E )C, transfer all storage capacities 
b* and the associated charging/ discharging profiles (7.*(t), 5*{t)), to the neighboring node j. In 
particular, consider the point (glit) , ^kit) , ^k{t), (^k{t),Pki{t), bk, k E M , k l,t E [T]j defined 
as follows. 





= 0, 7,w = 7;w+7;w, 




= 7^^), kEAr\{z,j}, 


kt) 


= 0, 6,it) = 6:it) + 6*{t), 


hit) 


= 6l{t), kE^\{z,j}, 


m 


= 9*{t) + -{^:{t)-6:{t)), 

Vij 


ek{t) 


= ei{t),kEM\{z}, 


h 


= 0, bj = b* + b*, 


bk 


= bl kE^\{z,j}, 


Pijit) 


=p:,{t)+i:{t)-6:{t), 


Pklit) 


= pUt), kr^l,{k,l)y^{z,J 



We do this successively for each i E IC to obtain a feasible point of U'^. Since the generation 



profiles remained invariant, the resulting point is optimal for H^. This completes the proof of 
Theorem [U 



B. Discussion 

First, we explore a few practical power networks, where Theorem \T\ applies, i.e., network 
topologies with generator buses that have single connections. In particular, consider the networks 
shown in Figure |3] The single generator single load case in Figure |3a] models topologies where 
generators and loads are geographically separated and are connected by a transmission line, e.g., 
see [|45l . This is common where the inputs for the generation technology (like coal or natural 
gas) are available far away from where the loads are located in a network. Figure |3b] is an 
example of a radial network, i.e., an acyclic graph. Most distribution networks conform to this 
topology, e.g., see ll38l . Il46l . Also, isolated transmission networks, e.g., the power network in 
Catalina island ffT3l are radial in nature. Theorem [T] also applies to generic mesh networks that 
have generator buses with single connections, e.g., the network in Figure |2l 

In all these examples. Theorem [U implies that for any available storage budget h, it is always 
optimal to place no storage on generator buses that have single connections. Problem P, in 
general, has multiple optimal solutions, but Theorem [T] proves the existence of an optimum with 
the property that b* = for all z in /C. 



Af\IC 




(a) (b) 
Fig. 3: Examples of power networks (a) Single generator single load system (b) A radial network. 



Our result is robust to changes in demand profiles, generation capacities, line flow capacities 



and admittances in the entire network, i.e., it remains optimal not to place any storage at buses 
in set /C under any changes to the above-mentioned parameters. The optimal sizing and the 
operation of the storage devices, however, might vary with changes in these parameters. To 
illustrate the efficacy of this robustness, consider the example in Figure |3al Suppose the line 
flow capacity is larger than the peak value of the demand profile, i.e., /12 > max^gj^^] d2{t). It 
can be checked that placing all the available storage at the generator bus is an optimal solution. 
If at a later time during the operation of the network, the demand increases such that the peak 
demand surpasses the line capacity, this placement of storage no longer remains optimal and 
requires new infrastructure for storage to be built on the demand side to avoid load shedding. 
If, however, we use the optimum as suggested by the problem II'^ and place all storage on the 
demand side from the beginning, then this placement not only can accommodate the change 
in the demand, but, it also, remains optimal under the available storage budget. To explore a 
different direction, suppose another generator is built to supply the load in Figure |3al Then, 
Theorem \T\ implies that the optimal placement still has no storage at bus 1 and thus is robust 
to such extensions of the network. In summary. Theorem \T\ suggests that the optimal storage 
placement problem has an interesting underlying structure and the solution to this problem is 
robust to changes in a wide class of system parameters and hence useful for network planners. 

The network planner uses Theorem [T] to solve the investment decision strategy and the optimal 
control policy for the generation and storage devices for a power network as follows. Let /C be 
the set of all generator buses with single connections. With the demand profiles and network 
parameters as input, solve the problem II'^ optimally. The optimal storage capacities bl,k e Af\lC 
define the investment decision strategy for sizing storage units at different buses. The optimal 
generation profiles gl{t), k E Mg define the economic dispatch of the various generators and the 
optimal charging/ discharging profiles 7^(t), 5l{t),k G A/" \ /C define the optimal control of the 
installed storage units. 

We end this discussion with a simple remark regarding the incorporation of renewables to 
the power network. Renewable generators with marginal cost of production can be treated as 
negative loads (and the corresponding buses as load buses). The result of Theorem 1 still applies 
to all conventional generator buses that have single links. However, buses with single connections 
with a renewable energy source differ from the corresponding ones with conventional sources 
and thus are not guaranteed to have zero storage in the optimal placement. 



C. On generators with multiple connections 

Generator buses with multiple connections may not always have zero storage capacity in the 
optimal allocation. In this section, we illustrate this fact through a simple example. Consider a 
3-node network as shown in Figure |4] All quantities are in per units. Let the cost of generation 
at node 1 be ci{gi) = gf. Let T = 4 and the demand profiles at nodes 2 and 3 be 

da = (9,10,0,10) and 4 = (0, 10, 10, 10). 

Also, suppose that the line capacities are /12 = /13 = 9.5 and the available storage budget is 
h = 5. Finally, assume no losses and ignore the ramp constraints in the charging and discharging 
processes, i.e. a = 1 and = es = 1. It can be checked that = 877 < irl^^ = 900.75. 



I generator bus 
# load bus 

3 

Fig. 4: A network with a generator that has multiple connections. 

IV. Results on specific network topologies 

In problems P and 11^, we solve for the optimal placement and control of storage in a power- 
network, given the demand profiles dk{t),t E [T], the storage budget h, the capacities of the 
generators g^.kE Mg and other network parameters such as the line flow limits fki, k I. In 
this section, we explore the behavior of the optimal cost of production as a function of these 
parameters. This provides valuable insights on various design issues, e.g., how much savings in 
terms of generation cost do we achieve by investing in an extra unit of storage. We explore such 
questions for specific network topologies. 

We make a few simplifying assumptions in this section. Let Ck{-),k E Mg be strictly convex 
and let a = 1 and = = 1. The proofs of the results are included in the Appendix. 

A. Single Generator Single Load Network 

Consider the single generator single load network shown in Figure |5l Generator at bus 1 is 
connected to a load (or demand) at bus 2 using a single line, i.e., /C = Mg = {1} and Md = {2}. 




9i{t)<9,® ® 

\Pl2{t)\ < fi2 



hi 



d2{t) 
r2{t) 



Fig. 5: Single generator single load network. Available storage budget h h>hi + 62- 



For this network, placing all the available storage resources at the load bus is always optimal. 
This is an immediate consequence of Theorem [Tl In this section, for any fixed demand profile 
d2{t),t E [T] of the load bus, we analyze the behavior of the optimal cost of production as 
a function of the generation capacity g^, the line flow capacity /12 and the available storage 
budget h; in particular, let the parameterized storage placement problem be P{gi, /12, h) and its 
optimal cost be /12, h). Similarly define, U^^^{'gi, /12, h) and ttI^^'Qi, /12, h). 

At the optimum of -P(^i, fi2, h), we have gi{t) < fi2,t E [T] from Lemma |2] Also, it satisfies 
dii'^) ^Qi^i^ [T]- Thus, to characterize the optimal point of P{g-^, fi2, h), it is equivalent to 
consider the constraint gi{t) < min {g-^^, /12} ,t E [T]. 

Proposition 1. For any h > 0, problem P{gi, /12, h) is feasible iff min {gi, fu} > /min, where 

fmin = max < max f '^'^=^ — ^LJ.^ j^iax [5i:z£i±l — | I ng) 

[i<i<rV t y ' i<ti<i2<T \^ t2-ti Jj 

Moreover, if min {^1, /12} > /mm, then p^{gi, /12, h) = p*(/mm, /mm, h). 

We interpret this result as follows. If either the line flow limit /12 < /min or the generation 
capacity g^ < /min, the load cannot be satisfied. Notice that fmin for /i > is no more than 
fmin for h = 0. Thus, storage can be used to reduce the cost of operation avoiding transmission 
upgrades and generation capacity expansion [,30 J . Interestingly, for /12 > /min and g^ > /min, the 
optimal cost of operation does not depend on the specific values of /12 and g-^^. From transmission 
or distribution planning perspective, investment in line and generation capacities over /min do 
not reduce the cost of operation. We provide an illustrative example at the end of this section. 

Next, we characterize the behavior of P{gi, /127 h) and its optimal cost p*{gi, /12, h) as a func- 



tion of h. For a given /12 and g^, the minimum required storage budget to serve the load depends 
on the demand profile d2{t),t E [T]. This may or may not be zero, depending on d2{t),t E [T], 
fi2 and g^. We calculate this minimum required storage budget, (say hmin) in Proposition [2l 
Also, it is easy to observe that as we allow larger storage budget, the generation cost does not 
reduce beyond a point, i.e., there exists hsat such that p^i^^, f 12, h) = p*{gi, fi2, hgat) for all 
h > hsat- We also calculate hgat in Proposition [21 First, we introduce some notation. Construct 
the sequence {Tm}m=o follows. Let tq = 0. Define iteratively: 

Tm= argmax , (17) 

T^-i+l<t<T y t - Tm~l J 

for 1 < m < M, where M is the smallest integer for which tm = T. Note that the sequence 
depends only on the demand profile d2{t),t E [T]. For any x G M, let [x]"^ := max(x,0). 

Proposition 2. Problem P{g^, f 12, h) satisfies: 

(a) If mm{gi, fi2} < maX(g[T] ^S^slMeI^^ fji^^ P(g-^, fi2, h) is infeasible for all h > 0. 

f )• P{9iJi2,h) is feasible iff h > h^in 

and fi2, h) is convex and non-increasing in h, where 



hmin — max 

0<t\<t2<T 



^ (d2(r) - min{^i,/i2}) 



Lr=ti + 1 

Furthermore, (^^ , /12 , h) is constant for all h > hsat, where 

t 7~m,_l 



(18) 



hsat = niax 

l<m<A/ 



max 

Tm-l+l<t<T„ 





(19) 



The condition min {g^, fi2} > maXf^^T] (^^^^^^~^^) implies that there is some /i > for 
which P(gi, fi2, h) is feasible. If this condition is violated, the problem remains infeasible no 
matter how large the storage budget h is. More the storage budget, lesser is the generation 
cost and hence p*igi, fi2,h) is decreasing in h. The convexity, however, implies that there is 
diminishing marginal returns on the investment on storage, i.e., the benefit of the first unit 
installed is more than that from the second unit. As a final note, observe that hsat is a function 
of only the demand profile and is independent of the generation and line flow capacities. 
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Fig. 6: Plots to illustrate propositions \T\ and l2l (a) Typical hourly load profile and optimal 
generation portfolio for line flow capacity /12 = 0.85, generation capacity gi = 1 and storage 
budget h = l (h) p^{g^ = 1, fu = 0.85, h). (c) p^{g^ = 1, /12, h = 1). 



Example: Now we explain propositions [I]and|2]with an example. All quantities are in per units. 
Consider an hourly load profile d2{t),t E [T] as shown in Figure |6al The optimal generation 
profile gl{t),t G [T] for P{gi = I, fu = 0.85, /i = 1) has been plotted in the same Figure. 
Notice that maxtg[r] gl{t) < /12 as stated in Lemma [2l 

Consider the plots in Figures |6b] and |6cl We plot p*{gi = 1, / = 0.85, h) for h in [0, 3] in 
Figure |6bl Notice that /12 < maxt^[T] d2{t), i.e., the problem is infeasible in the absence of 
storage. We calculate /imin = 0.226 and hgat = 2.598 from proposition |2l In Figure |6cl we 
plot p*{gi = 1,/12,/i = 1) for /12 in [0,2]. As in proposition [H the problem is infeasible for 
fi2 < /min = 0.683 and the optimal cost remains constant for /12 > /min- 



B. Star Network 

Consider a star network on n > 2 nodes as shown in Figure Ul J^g = {1} and A//? = 
{2, 3, ... , n}. For fixed demand profiles dk{t),t E [T],k E Md, line flow capacities fik, k E Md 
and capacity of the generator g^, let P{h) and n^^J^(/i) denote the storage placement problem 
and its restricted version as functions of the available storage budget h. Also, let p^{h) and 
'Kl^\h) be their optimal costs respectively. 




2 
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Fig. 7: A star network on n buses. 
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Fig. 8: P{h) and n^^^(/;,) for the simple 3-node star network in Figure S] 

In Section IIII-CI we showed that placing zero storage at the generator bus of a star network is 
not optimal, i.e., in general, p^{h) ^ vr* (/;,). In Figure [H we plot j>*(/i) and Tx^{h) for the 3-node 
star network shown in Figure |4] over a range of values of the total storage budget h. Observe 
that j5*(/i) < T\^^ ih) for some values of h but they coincide at: 

• Minimum value of h for which P{h) and Y\S^^{K) are feasible. 

• Large enough values of h. 

We formally state this for a general n-node star network in the following. Assume g = oo. 




■P(h) 



— ^^^y"^ — j fof <^ll k e Md- Then, P{h) and Il^^^{h) 

are feasible iff h > hmm, where 



max 

, 0<ti<t2<T 



t2 



E (4(r)-/u.) 



T=tl + 1 



(20) 



Moreover: 

(a) Pif{hrnin) ^* ij^min)' 

(b) There exists ho > /imm such that p^{h) = 7ri^^(/i) for all h > ho- 

V. Conclusions and Future work 

In this paper, we have derived analytic results on the optimal storage placement and sizing 
of storage units in the power grid. This provides valuable information for transmission and 
distribution system planners to find a suitable investment strategy in large-scale storage for the 
future smart grid. This work suggests that the storage placement problem has an interesting 
underlying structure and can be exploited to make sound planning decisions for sizing storage, 
that are robust to changes in a wide class of network parameters. As a result, it provides a 
framework to find the optimal storage capacity allocations in a network and in addition finds 
the optimal control policy for the generators and storage devices. 

A natural direction for future work is to study the same problem with performance metrics 
other than the production cost such as how storage can be used to defer transmission line 
upgrades or generation capacity expansion. This is a first step to "assess the technical and 
economic attributes" [|28l of this emerging smart grid technology. Another interesting direction 
is characterizing the optimal allocation of storage with stochastic demand and generation. In 
particular, we believe that storage placement on generators with multiple connections with 
realistic stochastic demand models would provide us with valuable insights. 
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Appendix 

Here, we present the proofs of the results presented in Section |IVl For the single generator 
single node and the star network, we drop the voltage angles 9k{t),k G Af,t e [T]. For any 
value of the power flow pik{t) from bus 1 to bus k, voltage angles 9k{t) can always be chosen 
to satisfy the power flow constraint in 

Furthermore, since a = 1, define rfc(t) := 7A:(t) — 5fe(t) be the power that flows into the storage 
device at node A; G A/" at time t G [T]. Notice that Vkit) can be positive or negative depending 
on whether power flows in or out of the storage device. Also, the storage level of the storage 
device at node G TV, at time t can be written as Skif) = Yl\=i'^'k{T)- 

A. Proofs for Single Generator Single Load Networks 

We drop subscripts from the variables d2{t), gi{t),t G [T], /i2,^i,ci(-) and the superscript 
from H^^^ {■) , ttI^^-) for ease of notation throughout this section. 

Proposition [B For any h > 0, problem P{g, f, h) is feasible iff min {g, f} > /min, where 

1 t6[T] V ^ / l<*i<fe<T \ t2-ti I \ 

Moreover, if min {g, /} > f^i,,, then p^{g, /, h) = p*(/min, /min, h). 

Proof: From Theorem [T] it suffices to show the claim for n(^, /, h) and 7r^(g, f, h). First, 
we show that if 11 /, h) is feasible, then min {g, f} > fmin- Fix any h>{) and let g{t),t G [T] 
be a feasible generation profile. Since X]t=i ^^{t) = S2{t) > 0, we have for any t G [T] 

t'G[Tr ^ ^ - t t - t 

Furthermore, for any I <ti < <T, the power extracted from the storage device between 



time instants ti and ^2 cannot exceed the total storage budget h and hence we have 

maxqlt > = > . (23) 

t'm t2-ti t2-ti - t2-ti 

Since g{t),t E [T] is feasible, g{t) < min{^, /} for all t E [T]. Hence, combining (l22l) and 
(|23]). we get 

min {g,f}> max g{t') > fmin- 

Next, we show that min {g, /} > /min is sufficient for n(^, /, h) to be feasible. Consider the 
optimal generation profile g*{t),t E [T] for the relaxed problem n(+oo, +00, /i). Suppose it 
satisfies 

max g*(t) < fmin- (24) 

*G[T] 

Then g*{t),t E [T] is also feasible and optimal for problem Il{g,f,h) for min {g,f} > fmin- 
Also, 7r*(^, /, h) = vr*(/min, /mill, h) for min {g, /} > /mm- It remains to show that (|24|) indeed 
holds. Consider the following notation. 



tie 



max{t E [T] \ g*(t) = max g*{T)}, 

{0 <t < tmax I g*{t) < g*{tmax)] ■ 



In the above definition g*{Q) := for convenience. If g*{tmax) = 0, then (|24|) clearly holds. 
Henceforth, assume g*{tmax) > 0. Then, g*{t),t E [T] satisfies: 



max q*(t) 
te[T] 



tmax 



tmax tless 

Et=wi [d{T)+r*{r)] 



tmax ti 



tmax tie 

Now, suppose the following holds: 



ess 



S2(tmax) ^2i.tless^ 



(25) 



S2{tr 



and s;(t;e.s) 



0, if tiess = 0, 
h, otherwise. 



(26) 



If (1261) holds, it follows from 

r=l d[T) 



^ ) if tless 0, 

maxg*{t) = ^ Y^^^^f d{r) - h 

' — , otherwise, 



tmax tless 

< f ■ 
— J mm- 

and hence (|24)) is satisfied. Next, we show that (|26|) indeed holds to complete the proof. First we 
prove that s^it^ax) = 0, i.e., the storage device at node 2 fully discharges at time t^ax- Suppose 
slitmax) > 0. As in Lemma [2l we construct a modified generation profile and storage control 
policy that is feasible and has an objective function value no greater than 7r*(+oo, +oo, h). But, 
the optimal generation profile g*{t),t E [T] is unique since the cost function c(-) is assumed to 
be strictly convex. Hence we derive a contradiction. By hypothesis, S2{tmax) > and hence the 
storage device at bus 2 discharges for some t > tmax- Let ti be the first such time instant and 
define 

Ai := mm {sKtmax) , g*{tmax) , 9* {tmax) - 9* {tl)] - 

Notice that Ai > 0. Consider the modified generation profile g{t) and control policy f2{t), that 
differ from g*{t) and r2(t) only at tmax and ti as follows: 

g(tmax) = g*{tmax) - Ai, ^(ti) = g* {ti) + Ai, 
flitmax) = rl{tmax) " Ai, r2(ti) = rl{ti) + Ai- 

Using Lemma [H we have 

c{g{tmax)) + c{g{ti)) < c{g*{tmax)) + c{g*{ti)). 

It can be checked that the modified profiles are feasible for n(+oo, +00, /i). The details are 
omitted for brevity. This is a contradiction and hence S2{tmax) = 0. 

Next, we characterize S2(t«ess)- If tiess = 0, then S2(tiess) = = 0. If tiess > 0, we prove 
that S2{tiess) = h, i.e., the storage device at node 2 is fully charged at time tiess- Suppose 
shitless) < h. As above, we construct a modified generation profile g{t) and storage control policy 
f2(t) that achieves an objective value no greater than 7r*(+oo, +00, h) to derive a contradiction. 



In particular, define 



A2 := mm{h - shitless) , g*(tless + 1) , g*(tiess + 1) - g*(tless)} > 0. 

Consider g{t) and f2{t), that differ from g*{t) and r2{t) only at tiess and tiess + 1 as follows: 

g{tless) = g*{tiess) + A2, g(tiess + 1) = g* {Uess + 1) " A2, 
f2{tless) = rliU^ss) + A2, f2{tiess + 1) = ^2(^«es. + 1) - ^2- 

As above, this defines a feasible point for n(+cxD, +00, h) and achieves an objective value strictly 
less than 7r*(+oo, +00, h). This is a contradiction and hence s*2{tiess) = h for tiess > 0. 



Proposition |2l Problem P{g, f, h) satisfies: 

(a) If mm {g, /} < maXig[T] ^^^^ -P(^7 f, h) is infeasible for all h > 0. 

(b) Suppose, min{^, /} > max(g[5"] ^ '^^'^^ j • T/zen, P{g,f,h) is feasible iff h > hj^i^ and 
p^:{'g, f, h) is convex and non-increasing in h, where 



^min — max 

0<ti<t2<T 



J2 {d{r)-mm{gj}) 



.T=tl+1 



(27) 



Furthermore, p^{g, /, h) is constant for all h > hgat, where 

t ^m — 1 



hsat = max 

l<m<M 



max 

Tm-l+l<t<T„ 



\T = Tm_l + l 



(28) 



Proof: From Theorem [T] it suffices to prove the claim for Il{g, f, h) and 7r*((7, /, h). 

(a) To the contrary of the statement of the Proposition suppose that min {g, /} < max^gfy] {^^^^^-^ 
and n(^, /, /;,) is feasible for some /i > 0. Then, it follows directly from Proposition [Tlthat 
min{^, /} > frnin > maxte[T] contradicting our hypothesis. 

(b) First we show that if 11 /, h) is feasible then h > hmin- Suppose 11 /, h) is feasi- 
ble. Then, for all < ti < ^2 < 7" Proposition \T\ implies that min {g, f} > fmin > 
(Er=ti+i d{r) - h) /{t2 - ti). Rearranging this we get h > Y!r=H+i (^(^) " ^in {g, /}) . 



Also, h >0 and hence: 

t2 



h > max 

0<ti<t2<T 



^ {d{T)-mm{9j}) 



+ 

= h 



mm- 



LT=tl+l 

Now we show that h > hmm is sufficient for n(^, /, h) to be feasible. The relation h > hmin 
can be equivalently written as follows: 

V*^ , d(T) - h 

min {g, /} > ^"-*^+^ ^ ^ , for all 0<ti<t2<T. (29) 

to — U 



Also, by hypothesis, we have 



min / } > max ( ^r=id{r) , _ ^^^^ 



te[T] \ t 

Combining (|29l) and (|301 ). we get min /} > /,™n- Then, Proposition [T] implies that 
n(^, /, h) is feasible. Convexity and non-decreasing nature of p*{g, f, h) as a function of h 
follows from linear parametric optimization theory [ )39l . 

Finally, we prove that p*{g,f,h) is constant for all h > hsat, where hgat is as defined 
in (1281) . The proof idea here is as follows. We construct the optimal generation profile 
g*it),t E [T] for the problem n(+oo, +00, +c)o) and show that it is feasible and hence 
optimal for the problem n(^, /, +00) provided min{^, /} > ma,xti=[T] hoids. 
Problem n(+oo, +00, +00) can be re-written as follows. 

T 

minimize y^c-i(q(t)) 

t 

subject to g{t)>0, ^(^^(r) - ci(r)) > 0, t e [T], (31a) 

T=l 

T T 

'£g{r) = J2d{r). (31b) 

T=l T=l 

Let the Lagrange multipliers in equations (|31al) - (|31bl ) be A(t), £{t), t E [T] and u, re- 
spectively. It can be checked that the following primal-dual pair satisfies the Karush-Kuhn- 
Tucker conditions and hence is optimal for the convex program n(+oo, +00, +00) and its 



Lagrangian dual [|39l . We omit the details for brevity. 

y^Tm d(T) 

9*{t) = ^"="-^+^ , t = r„„i + l,...,r^ and m = l,2,...,M, 

\c'{g*{Tj)-c'{g*{T^ + l)), ift = Tm, m = l,2,...,M-l 

r(t) = I , te [T], 

I 0, otherwise 

A*(t) = 0, t G [T], and z/* = -c'(^*(T)). 

The above profile g*{t),t G [T] of n(+oo, +00, +00) satisfies: 

max q'*(t) = max — ^^^^^ < min |q, f | , 

and hence is feasible and optimal for n(^, /, +00). Note that Xlr=r,„_i+i i9*i'^) ~ '^i'^)) = 
for all 1 < m < M. Thus, for r^-i < t < Tm, we have 



E"'" ^1 d(T) 



"^m ''"m— 1 I -, 

T=Tm-l + l 



Maximizing the above relation over all t G [T] we get maxjg[T] S2(t) = h^at- Therefore, 
g*{t),t G [T] is feasible and optimal for n(^, /, h) provided that h > hgat- 



B. Proofs for Star Networks 

— j for all k G Md- Then, P{h) and Tl^^^Qi) 

are feasible iff h> hmin, where 

t2 



hmin — 7 , max 

^ 0<ti<t2<T 



iMr)-fik) 



.T=tl + 1 



(33) 



Moreover: 

(a) Pif{hmin) ^* {hmin)t 

(b) There exists ho > hmin such that p*{h) = nl^^h) for all h > ho- 



Proof: First we show that h > hmin is necessary for P{h) to be feasible. Consider any 
feasible solution of P{h). For any k G Md and < ti < t2 < 7", we have Y^T=ti+i ^ki^) — ~^ki 
since the power extracted from the storage device at node k cannot exceed the corresponding 
storage capacity hk. Also, for any k G Md the power flow on the line joining buses 1 and 
k satisfies pik{t) = dk{t) + rfc(t) < /i^ for all t G [T]. Combining the above relations and 
rearranging, we get bk > Xlr=ti+i {dkij) — fit)- Also for k G A/d, bk > and hence 

+ 



bk > max 

0<ti<t2<T 



.r=ti + l 



(34) 



Thus we get h > J^keAf^ ^k ^ ^mm- If Il^^^{h) is feasible, then P{h) is also feasible and hence 
h > hrnin IS ncccssary for both problems to be feasible. Now we prove that it is also sufficient. 
In particular, we show that for h = hmin, Ii^^^{h) is feasible. For convenience, define 

- t2 ] + 

k G Afo- (35) 



max 

o<ii<t2<r 



-=ti+i 
Then hmin 

= T,keMD Rearranging (|35]), we get 



fik > max . (36) 

0<ti<t2<T \ t2-U 



/i,>max(£^^). (37) 



Also, by hypothesis, we have 



Combining equations (l36l) and (l37l) . we have 

> max l^^max^^ ^ j , max ( , ) } . (38) 

For each k G Md, consider a single generator single load system as follows. Let the demand 
profile be dk{t), the capacity of the transmission line be fik and the total available storage 
budget be hk- For this system, the right hand side in (|38T ) coincides with the definition of fmin 
in (fT6l) . From Proposition [T] it follows that there is a feasible generation profile (say g'^'^\t)) and 
a storage control policy rk{t) that define a feasible flow over this single generator single load 



system and meet the demand. Now, for the star network, construct the generation profile gi (t) 

9i{t) = E 

and operate the storage units at each node k G Md with the control policy rk{t) defined above. 
Also, ri(t) = for all t E [T]. It can be checked that this defines a feasible point for Il^^^{hmin)- 
Next, we prove that p*{hmin) = T^l^Hhmin)- Let bl,k E Af be optimal storage capacities for 
problem P{hmin)- Then the optimal storage capacities satisfy the following relations: 

^hl> hmin, and 6* + ^ bl < hmin- 

where the first one follows from (|34T l and the second one follows from the constraint on the 
total available storage capacities. Rearranging the above equations, we get bl = and hence 
P*{hmin) = 7r=l^^(/imm)- TWs complctcs the proof of part (a). 
To prove part (b) of Proposition |3l we start by showing that 

p.(oo) = 7r«(oo). (39) 

Assume P(oo) is feasible. For /i = oo, we drop the variables bk,k E N, and consider the 
problems P(oo) and n''^^J^(oo) over the variables gi{t), rk{t),k E Af. The variables pik{t) and 
Sk{t) are defined accordingly for all k E M . We argue that the optimal points of P(oo) lie in 
a bounded set. Note that \pik{t)\ = \dk{t) + rk{t)\ < fik and thus the control policies rk{t) are 
bounded for all k E Md- Also, the cost function ci(-) is convex and hence its sub-level sets 
||39l are bounded. From the above arguments and the power-balance at bus 1, the optimal policy 
ri{t) is also bounded. Thus, the set of optimal solutions of P(oo) is a bounded set. Furthermore, 
this set is also closed since the objective function and the constraints are continuous functions. 
As in the proof of Lemma |2l associate the function Xlt6[r] ^^^^ every point in the set of 

optimal solutions of P(oo). This is a continuous function on a compact set and hence attains a 
minimum. Consider the optimum of P(oo) where this minimum is attained. We prove (|39] ) by 
showing that r*(t) = for all t E [T] at this optimum. 

Assume to the contrary, that r*(t) is non-zero for some t E [T]. Define 

to := {t G [T] I ri(to) > 0} and ti := min {t G [to + 1, T] | r*(t) < 0} . 



Also, define A := min {r]'(to), — ^i(^i)} and notice that A > 0. 

Case 1: glito) > gl{ti) + A: Construct the modified generation and charging/ discharging 
profiles gi{t),fi{t) that differ from gl{t),r1{t) only at to and ti as follows: 

gi{to) = g*i{to) - A^, ^i(ti) = glih) + A^, 
ri(to) = rlito) - Ag, hih) = r*(ti) + Ag, 

where Ag := min {A, gi{to)} > 0. As in the proof of Lemma|2l this is feasible for P(oo). Also, 
by Lemma [T] 

ci (^i(to)) + ci (giih)) < ci (glito)) + ci (gHh)) . 
The details are omitted for brevity. This feasible point satisfies 

|ri(to)| + |ri(ti)| = r*(to) - A - rUh) - A < \rl{to)\ + \rl{h)\, (40) 

and hence defines an optimal point of P(oo) with a strictly lower value of the function XltG[T] I '"i ('^) I • 
This is a contradiction. 

Case 2: gKto) < gl{ti) + A: As above we construct modified storage control policies ffc(t) 
for all k E J\f, keeping the generation profile constant to define an optimal point of P(oo) with 
a lower value of J2te[T] derive a contradiction. 

Let the modified control policy at bus 1 be as follows: 

h{to) = rlito) - A, fi(ti) = rt(ti) + A. 
Instead, we distribute this to storage devices at A; G Md, as follows: 

fk(to) = rl(to) + ^Jk, fkih) = rliti) - i^k, k e Md, 

for some V'/fc ^ 0, A; G Md and J2k£j\fo '^^ ~ ensure feasibility of the modified profiles 

it suffices to check that the line flow constraints are satisfied at to and ti. In other words, we 
show that there exists ipk^k E Md such that for all k G Md, 

i^k > 0, Pu-(^o) + i^k< /ifc, plki^i) -i^k> -fik, ^ v^fc = A. 



Equivalently, we prove that 



^ mill {fik - p*ik{to), fik + P*ik{ti)} > A. 

Recall that plf,{to) and Pifc(ti) are feasible for P{oo). Thus Pifc(to) < fik and plk{ti) > —fik- 
Also, gl{t) - rl{t) = Y.keMn'PikiP) at t = to and t = h. Thus, we have 



J2^^{fik-p*ikito),fik+plkiti)}> Yl (Plkih) - PIM) 




>-A <-A >A 



> A, 

where the last inequality follows from the hypothesis gl{to) < gl{ti) + A. The modified profiles 
satisfy |ri(to)| + |^i(ti)| < ki(^o)| + ki(^i)l as in (|40|) . As argued above this is a contradiction 
and hence (l39l) holds. 

For P(oo), sl{t),k e Af,t e [T] is finite. Define ho := XlfeGATo "^^^telT] Then, note 

that {gl{t),rl{t),t e [T] k e TV) are also feasible for U^^^h) and P{h) for all h > ho- This 
completes the proof of Proposition |3l ■ 



